Role of domain wall fluctuations in non-Fermi liquid behavior of metamagnets 
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We study resistivity temperature dependence of a tiiree dimensional metamagnet near the meta- 
magnet phase transition point in the case when magnetic structure tends to split into regions with 
high and low magnetization. We show that in the case of weak pinning the spin relaxation time of 
domain wall is much larger than that of the volume spin fluctuations. This opens a temperature 
range where resistivity temperature dependence is determined by scattering of conducting electrons 
by the domain wall fluctuations. We show that it leads to quasi-linear low temperature dependence 
of resistivity. 

PACS numbers: 



I. INTRODUCTION 

Understanding the deviation from the Fermi liquid be- 
havior at phase transition critical points is a current 
research question. The most important are deviations 
from quadratic temperature dependence of resistivity 
near quantum critical points. Theory proposes that ori- 
gin of non Fermi liquid behavior is scattering of conduct- 
ing electrons by bosonic critical soft modes. In case of 
ferromagnetic transition these modes are collective spin 
excitations virhose relaxation time diverges near the tran- 
sition point. In nearly magnetic metals at temperatures 
larger than inverse spin relaxation time resistivity due 
to electron scattering by spin fluctuations, strongly devi- 
ates from Fermi liquid quadratic dependence. It becomes 
linear or even saturates ji^— In strongly Stoner enhanced 
paramagnetic metals deviation might starts at very low 
temperature. 

It was experimentally shown that in the region near 
quantum critical point, magnetic state of many systems 
experience first order phase transition with separation of 
different phases.-"^^ In this case the temperature depen- 
dence of resistivity is characterized by non-Fermi liquid 
resistivity exponents, which depend on proximity of the 
system to the quantum critical point. And it was shown 
that dynamics of magnetic state in this case is not criti- 
cal. 

Nonuniform magnetic state can be considered as a dis- 
tinctive phase itself ji^ii^ For example in metamagnetic 
transition of thin magnetic films magneto-dipole interac- 
tion can give rise to formation of magnetic domains. In 
another example domains can be formed in spatially ran- 
dom variation of effective magnetic field. This model was 
proposed to explain unusual magnetotransport effects of 
ferromagnetic alloys 

The phase separation in first order magnetic phase 
transition might exists in considerable range of tempera- 
ture, pressure and other parameters. In broadened phase 
transition the ratio of phases volume, structure and area 
of phase boundary depend on temperature. In broad- 
ened phase transition such factors as large difference in 
phases resistivity, or strong electron scattering by the 



phase boundary can strongly influence temperature de- 
pendence of electronic kinetic properties of material. 

In this paper we consider the temperature dependence 
of resistivity due to scattering of conducting electrons 
by domain walls fluctuations in broadened metamag- 
netic first order phase transition. Domain walls corre- 
spond to crossover from states with low magnetization to 
high magnetization. Theoretically domain walls between 
phases appear in regions of space where two minima of 
free energy have equal value. The contribution to resis- 
tivity will only be significant when the wall fluctuations 
relaxation time is comparable or larger than that of the 
volume magnetic fluctuations. 

The paper is organized as following. In section II we 
give a description of the system of coupled conducting 
electrons and itinerant electrons responsible for metam- 
agnetic state. We propose a model where the magnetic 
state splits in to domain walls. In section III we obtain 
solution for the domain wall profile and consider the fluc- 
tuations around it. We then obtain dynamical magnetic 
susceptibility of the domain wall fluctuations. In section 
IV the temperature dependence of the resistivity caused 
by the scattering of conducting electrons by the domain 
wall fluctuations is considered. We find transition from 
quadratic to linear temperature dependence with increas- 
ing temperature. While this temperature dependence co- 
incides with that due the volume spin- fluctuations.^, it 
starts at much lower temperature. In section V we dis- 
cuss the interaction correction^* at small temperatures. 
It is known that conducting electrons scattering by impu- 
rities and spin fluctuations results in important contribu- 
tions to the temperature dependence of conductivity.^^" — 
We obtain that for domain wall fluctuations the interac- 
tion correction to conductivity has temperature depen- 
dence predicted for the 2D volume spin fluctuations <^ 

II. DESCRIPTION OF THE MODEL 

We consider a three dimensional metallic metamagnet 
where s— electrons are considered to be conducting, and 
d— electrons responsible for the magnetic state. Coupling 
of conducting electrons with bosonic modes is obtained 



by integrating out d— electrons ^23. Action describing con- 
ducting electrons in random impurity potential V(r) is 
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where "0 describes electrons with spin a and mass m, 
Matsubara frequency f2„ = i:T(2n + 1) (T is tempera- 
ture), with Fermi level /i, and we have set T% — 1. We 
assume impurity potential to satisfy < V{r) >= 0, and 
< V{r)V{r') >= l/{2Tn^T)6{r - r'), where v - density of 
electrons per one spin, r - is the electron mean free time. 

The free energy of the magnetic state is described by 
two energetically unequal minima. Application of mag- 
netic field brings the system in to a state at which these 
two minima have the same energy, this point is called 
metamagnetic phase transition point and the value of 
the magnetic field at which it occurs is called metamag- 
netic value of magnetic field. At the metamagnetic phase 
transition point we will be describing the free energy of 
the metamagnetic state by two parabolas with minima 
at TO(r,T) = ±mo, corresponding to high and low mag- 
netization states. In this case the domain walls originate 
from the spatial deviation of the magnetic field from its 
metamagnetic value, and we introduce such deviation as 
h{r) (see figure 1). Assuming strong magnetic field we 
only consider longitudinal component of magnetization 
density m(r) (in units oi gfiB = 1) in the action. Having 
above assumptions in hand we write the action describing 
the magnetization density m(r, r) in the form 

S.ma,H^Jo^'^dTjdv[^jWm{r,T)f 
+ 2S_(|m(r,T)|-mo)2-/i(r)m(r,T)) (2) 

here coefficient K~^^^ is of the order of electrons Fermi 
wave length, a^^ is a Stoner enhancement factor, and xo 
is conducting electron spin susceptibility which we have 
assumed to be equal for both high and low magnetiza- 
tion phases. This assumption greatly simplifies further 
calculations, and does not affect main conclusions about 
resistivity temperature dependence. We assume the spa- 
tial distribution of the h{r) to be of the Gaussian form 
and give the averaging procedure in the appendix. 
The last term in ^ is the Landau damping described 

by 
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dr'dr to (r, cl;„) r(r - r', a;„)m (r', a;„) 



(3) 

here w„ is the bosonic Matsubara frequency. At this 
point it is necessary to distinguish two cases of the form 
of the damping. When the scattering of conducting elec- 
trons by impurities is ballistic the Fourier image of F is 
given by F(Q,cj„) — ^^"g , which is valid for large mo- 
menta vfQ > |w„| where vp is Fermi velocity, and 7 
is a damping constant. ^'^^ In case of diffusive scatter- 
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FIG. 1: Schematic phase diagram of the broadened first order 
phase transition. The solid curve shows putative transition, 
while area enclosed by dashed curves represents the phase 
separated state. 



which is valid for 



ing the damping is F(Q,i:j„) — -pg^ 

DQ^ > \Ldn\ where D is the diffusion constant;^ 

Finally, action describing coupling of conducting elec- 
trons with magnetization is 



Sint = ^ y" dr J dr s{r,T)m (r, r) , 



(4) 



here s (r, r) is operator of spin density of conducting elec- 
trons along the longitudinal component of magnetization, 
G is a phenomcnological coupling constant. 



III. DOMAIN WALL FLUCTUATIONS 

In this section we give a solution to the mean field 
equation for the domain wall and discuss fluctuations 
around this solution. Let a; be a coordinate normal to 
the domain wall and h{xo{p)) = 0, so that at x < xq 
it is state with high magnetization, and at x > xq it is 
a low magnetization state. Here p is a 2D coordinate 
along the two dimensional domain wall. Equation for 
magnetization is obtained by varying the action ^ and 
we get 



K -—rm + a (m — mosignirn)) — xoh(x) . (5) 
dx'' 



When h{r) is slowly varying function on a scale of K/a 
the domain wall can be approximated as flat. With this 
assumption the solution of (O describing the domain wall 
is 

Xoh{x) 



a 



TO = -TOo(l-e"V^I=^-^o(p)l)gj^„(2,_2;o(p)) + 

( V 

(6) 

Let us consider fluctuations near this solution. Taking 
second derivative of free energy we obtain equation for 
eigenfunctions, which describe fluctuations 

{-KV^ + a)Sm{r) 
2amo S{x-xo{p))6m{r) ^e6m{r). (7) 



3 



Deriving this equation we have used the equahty 

6{x - xo{p)) 



6{m{x)) 



|;f(m(a;))| 

Idx^ y >>\ x=xo 



(8) 



S type potential in equation ([7]) is related to non ana- 
lytical dependence of free energy on magnetization. The 
equation ([7]) has only one bounded solution, thus strongly 
simplifying consideration of fluctuations. At slowly vary- 
ing xq{p) and ■^m{xo{p)) one can search for the solution 
of equation l^i in the form of a plane wave in p 



5m{r) = vl/Q(a; _ a;o(p))e 



where 
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and assuming e = KQ^ + eo wc find 



eo = a 
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where /3 — {a — eq) / K . At slowly varying h{v) we find 
that eo << a and we approximate (3 — ^Jaj K. In this 
case 



2X0 

mo/3 



dx 



h{x) 



(12) 



Dynamics of spin fiuctuations is governed by the Lan- 
dau damping ([3]) which for excitation described by 5*0 (a^) 
translates to 



7(0) - - 



dxdx''^ o{x)'^oix') 
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where as an example we used ballistic case. At small 
momenta {3 > Q we have for the ballistic case 



7(Q) 
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Same procedure for the diffusive case gives 
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7(0) = 
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The dynamical susceptibility of one domain wall fluctu- 
ations is represented in the form 



X(r,r','^) 
where 



i!Q_e-*Q(p-p') 



^o{x) *o (x')x(Q,^), 
(16) 



x(Q,^) = 



Xo 



eo + KQ^ + iuj (Q) ■ 



(17) 



We consider a locally fiat domain boundary. This can 
be justified if scale along wall Ln ^ 



much smaller than domain wall curvature. Typical scat- 
tering of conducting electrons by fluctuations of the do- 
main wall is characterized by the momentum of the order 
of /3. Which is much larger than the momentum of the 
dynamical susceptibility of fiuctuations. This allows to 
separately average over the domain wall direction and 
position, introducing concentration of wall nw- The av- 
erage over random eo is approximated by substituting the 
average value of eo in to the susceptibility. The details 
of the described above averaging steps are given in the 
appendix to the paper. 

From (fT7|) we see that relaxation time of domain wall 
fiuctuations is proportional to e^^ which is much larger 
than that of the volume fluctuations since eo << ct. At 
small eo contribution of the domain wall fluctuations to 
the total susceptibility of the system is approximately 



6x0 ^nw dxdx'i'o (x) *o {x') x(0, 0) 



, (18) 



which can be of the same order as xo depending on the 
parameters. 



IV. RESISTIVITY TEMPERATURE 
DEPENDENCE 

Let us consider contribution to the resistivity due to 
domain wall scattering of conducting electrons. There 
are three temperature dependent contributions to the 
conductivity. They are due to scattering of conducting 
electrons by fiuctuations of domain wall, due to domains 
walls shape change when temperature increases, and due 
to variation of concentration of walls nw ■ First two con- 
tributions are considered in this section and have a com- 
mon nature. Third contribution depends on position on 
the phase diagram. We discuss its contribution in the 
conclusions of the paper. 

Resistivity due to electron scattering by spin fluctua- 
tion is obtained in second order perturbation theory in 
interaction Q and is expressed through imaginary part 
of averaged susceptibility as^i^ 

2pF 3 ^ 

p(T) = Rn— f —^-7— f ^ ImY(q,[j), 

' T J p% J sinh2(a;/2r) ^ ^' 
—00 

(19) 

here 
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where u is density of states of conducting electrons per 
spin, pp, n, m are Fermi momentum, density, and mass 
of conduction electrons respectively. 

In addition the fiuctuations give temperature depen- 
dent contribution to the average magnetization of wall. 
Fluctuation part of magnetization is determined by 
derivative dm{r) — — of fluctuation part of the free 
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energy Afl = jT^^c^^.q ln(eQ + j{Q)uJn)- This leads 
to a change of the domain wall profile and gives addi- 
tional temperature dependence of resistivity proportional 
to G^m{Q)Sm{Q,T). Sum of both contributions to the 
resistivity is given by 



rt — 1 



oo , 

f duj ( ^ . , :! 
n ^ 



("/2T) 



2coth(^) + 2 



/^Imx(Q, 



(21) 



Term, proportional to (— coth(^) + 1), is due to 
(5m(Q,r). 

Our calculations show that despite of the difference 
in damping in ballistic ([T4|) and diffusive (jisp regimes 
the temperature dependence of resistivity has the same 
form for both of the cases. It is quadratic at tempera- 
tures lower than the inverse relaxation time and linear at 
larger temperature. We get that in the ballistic scattering 
regime at temperatures T < Tq resistivity has quadratic 
temperature dependence 
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(22) 



and at temperatures T > Tq the dependence is linear 



R = 



2TTl3RonwXo T 



(23) 



here Tq — /3wFeo/(5 ln(/3^iir/eo)). In case of diffusive 
scattering regime we get that the resistivity dependence 
of the temperature has the same form with Tq modified 
as To = Dl3ei^\/ K / tQ. The transition from quadratic 
to linear T dependence corresponds to transition from 
quantum to thermal fluctuations of domain walli^ 

We would like to notice that the same temperature 
dependence holds for scattering by volume spin fluctu- 
ations, except for the difference of effective To in bal- 
listic and diffusive cases. ^'''^^ Let us compare resistivity 
temperature dependence (|22p and that Qvoi{T) due to 
volume spin fluctuations. In the considered temperature 
range the volume contribution is quadratic and Qyoi (T) ^ 
RoXoTyipFVFV^)-'''^ Ratio R{T)/g,oiiT) ^ ^ can 
be of order of one. 



V. INTERACTION CORRECTION 

At low temperatures important resistivity temperature 
dependence is related to weak localization and electron 
interaction corrections.— Here we are going to discuss 
contribution to the conductivity originating from the in- 
terplay between electron inelastic scattering by domain 
wall fluctuations and elastic scattering by impurities. 
The triplet channel contribution to the conductivity after 
disorder averaging is given by'^^' ™ 



6a = 2ne'vlr,yG' 1 [£ coth ^)] 



(24) 



here we also averaged over the positions of the do- 
main walls after which the equation above became 
isotropic. In the ballistic regime Tr > 1 one can 
use the following approximation i?(q, cj) « 2/{vpqY 
which is valid for vpq > and the expression 

for the domain wall profile fluctuations propa- 
gator at small momenta q < (3 is x{''\i^) = 

Wxo//3 (eo + KQ'' + / {^vfP)) In {a/KQ-')y\ 
We only consider temperature dependent terms of the 
correction. 

We find that at temperatures T > TrvFl3eo/{2'-f) ~ Tq 
and T < vp^/^a/K the correction to the conductivity is 
logarithmic in temperature 



Sa 



f Xq 



\2Tr^j 



e^vG^n^Y ) a/ eo/KvpTln 



(25) 

At higher temperatures T > vp^JT^fK and T < 
the correction to the conductivity is linear 



(26) 



At higher temperatures the correction decays as l/T. 

Now let us discuss the diffusive regime at Tr < 
1. In this case B(q,a;) = I (D$\i^ f 

X(q, t^) = 4nw'Xo//3 (eo + if 0' + j27u;/(^/30))7' ap- 
proximated at small momenta q < /3. Most singular 
contribution arises from ^/\uj\/D < q < min(/3,l/^), 
and q > y^eo/K. Where £ is the electron's mean 
free path. Calculations show that at temperatures 
p = min(r2i, ri2, l/r) > T, where fli = ^ 
^/2f3'^DK/j the contribution to the conductivity is 
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27/2 



^e^i^G'^nw j In (^^^ In 



pT 



(27) 



We would like to point out that we have not considered 
all of the possible regimes of parameters, focusing on the 
most interesting cases. Expressions (^5]) . (pSj) . and (P7)) 
have a 2D temperature behavior while obtained for the 
3D electron system.^ The resistivity correction is ob- 
tained by 6p = 1/(7 £i(l — Sa/ao), where an is the Drude 
conductivity. 



VI. CONCLUSION 

We considered the temperature dependence of resis- 
tivity due to electron scattering by spin fluctuations of 
domain walls. We showed that in case of weak pin- 
ning a >> eo the relaxation of domain wall fluctuations 
is much slower than relaxation of volume spin density. 
Therefore, classical regime of fluctuations characterized 
by linear temperature dependence starts at much lower 
temperature than for the volume fluctuations. In consid- 
ered temperature range contribution of volume fluctua- 
tions to resistivity is always ~ T^. 
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It is reasonable to assume, that additional resistivity 
in considered system is proportional to n^. Depending 
on relation between average magnetic field and metam- 
agnetic hm (T) the concentration n^, can change with the 
temperature. In considered temperature range the con- 
centration changes with the temperature as nwiT) = 
nvF(0)(l + a{TTs)'^), where \a\ ~ 1 and the sign of 
a depends on the average magnetic field. The ts ~ 
1/ (appvp) is relaxation time of volume fluctuations. 
Parameter Tts determines renormalization of free en- 
ergy due to volume spin fluctuations."^ For example in 
schematic phase diagram shown in fig. 1 the nyy de- 
creases in the region h < when h„i{T) increases with 
temperature. 

In case when a < there can be a cancelation of 
quadratic temperature dependence of volume contribu- 
tion to resistivity by the contribution of domain walls 
nw{T). In case of the cancelation, obtained results (P^. 
(PSj) . and ([25]) will be dominant. And the total resistiv- 
ity will have quasi-linear non-Fermi liquid temperature 
dependence. 
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Appendix A: Averaging over random surfaces 

Here we show the averaging procedure over the direc- 
tion of domain walls and random cq. We assume that 
magnetic field is described by Gaussian distribution such 
that 

(Aft(r)A/i(r')) = Aexp(-C(r - r')^) (Al) 

and (A/i(r)) = 0, where Ah{r) = h{r) — Hq. When pa- 
rameter C is small the magnetic field becomes a slowly 
varying random function h(r). Typical phase boundary 
in such a random field is smooth. Position of domain wall 
is defined by the following equation 

h{r) =hQ + Ah{r) = (A2) 

Note, that under this definition we neglect cases when 
in small closed regions there is no solution of mean field 
equation ^ for favorite phase, or energy associated with 
it is too high. 

Fourier transform of any quantity V{r) which is 
nonzero near the domain wall surface and slowly vary- 
ing along it, is calculated as 



V{q) = J dr exp(iqr)l/(r) = J dS exp(iqrs)F(qn(rs), rg) ^ J dr l/(qn(r), r) exp(zqr)|^/i(r)|^(/i(r)), (A3) 

here is a point on the surface, and n(r5) is a vector normal to the surface at the point rs defined as n(r) = 
^^^^/|^/i(r)|. In (jA3l) V^(qn(r),r) is a one dimensional Fourier transform in the direction normal to the surface. 

We now can average the dynamical susceptibility of the number of domain wall. For one domain wall the sus- 
ceptibility is given by expression (jl6p . Scale along the surface of the domain wall that we are interested in is 

Lii . 1^ — \= here 5h ~ ^/~A ~ and the en is estimated as 

II V '^0 PF\/C0' XO 

^o-^^lf^l-v/^yC. (A4) 
mo V a/if dx 

Therefore Ly ~ j^VP/'^VC is smaller than the radius of surface curvature i||-\/C ^ V V^/f^ << 1; and under 
averaging procedure we can use the expression (|16l) . Then the average susceptibility is 



X(q,c.) ^ Jd^'^J rfbn(R,b)exp(zqR)b2«'o(qn)*o(-qn) J exp(iQR)x(Q, (A5) 

where ^o{q) = 2/33/2/(^2 ^ ^2) 

is a Fourier transform of ^o(a;) defined by PH)) . and 
n(ri - r2, b) = {6{h{r,))Sihir,))Sih - \{^h{v,) + ^Kv^- (A6) 

I 



At i? < < the b is a normal to the surface and there- 



fore is also normal to R . Therefore in (lASp we can sep- 
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arately average over the direction and value of b, consid- 
ering n(R, b) only in the limit of Ry/C, « 1. In case of 
Gaussian distribution n(R, b) is calculated analytically 

as 



where 



and 



nw 



n(R, b) = ^^F(b^)P(b|| ), (A7) 



P(b|| 



exp(- 



3b^ 



(A9) 



here bj^ and by are perpendicular and parallel compo- 
nents to R consequently. The quantity nw has a mean- 



ing of the domain wall concentration and is defined as 



nw 



{Jdr6ihir))\£h{r)\) 



V 



-^exp(-M) (AlO) 

TT A 



where V is the volume of the system. The factor b 



under the integral of expression (jASP is estimated as b ~ 
Eq. So averaging of quantities like susceptibility does not 
diverge at small eg and can be approximated by inserting 
average eq. Finally, the averaged susceptibility over the 
domain wall positions is given by the next expression: 



^ J d2n*o(qn)^'o(-qn) 



nwXo 



(All) 

where ep is approximated by its average value. At small 
momenta q < /3 Fourier transform \I'o(qn) is approxi- 
mated as *o = 
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